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ABSTRACT 

We present results of a Chandra survey of the ultra-luminous X-ray sources (ULX) 
in 13 normal galaxies, in which we combine source detection with X-ray flux mea- 
surement. 22 ULX were detected, i.e. with L x > 1 x 10 39 erg s _1 (Xio)> and 39 other 
sources were detected with L x > 5 x 10 38 erg s^ 1 (L 5 ). We also use radial intensity 
profiles to remove extended sources from the sample. The majority of sources are not 
extended, which for a typical distance constrains the emission region size to less than 
50 pc. X-ray colour-colour diagrams and spectral fitting results were examined for 
indicators of the ULX nature. In the case of the brighter sources, spectral fitting gen- 
erally requires two-component models. In only a few cases do colour-colour diagrams 
or spectral fitting provide evidence of black hole nature. We find no evidence of a 
correlation with stellar mass, however there is a strong correlation with star formation 
as indicated by the 60 /xm flux as found in previous studies. 

Key words: accretion: accretion discs - binaries: close - black hole physics - X-rays: 
galaxies - X-rays: stars 



1 INTRODUCTION by a multi-colour disc blackbody model suggesting black 

hole nature. Chandra observations of galaxies have revealed 

TT1 , , . v / TTT -v\ j. ■ ■ a large number of previously unknown ULX, the luminosi- 

Ultra-luminous X-ray sources (ULX) are intriguing, appar- 6 * J 39 -i ^ 

ties ol which extend from 1 — 10 x 10 erg s (Fabbiano, 



Zezas & Murray 2001; Blanton et al. 2001), suggesting that 



ently point-like sources in external galaxies which are nor- 
mally distinguished from bright central objects. Bright ob- 

. , , ,. , . , ... , they may be quite common. With the advent of Chandra, it 

jects may be divided into those that exceed the Eddmg- J , 

c , . „ n . .,„38 becomes possible to study ULX in more detail, 

ton limit LEdd for a 1.4 Mq neutron star of ~ 2-4x10 J 

erg s _1 , and those with luminosity greater than 10 39 erg The discovery of ULX naturally led to the proposal that 

s , which is the normal definition of an ULX. Observations these may be a single type of object. Explanations of this 

of bright spiral galaxies with Einstein revealed substantial type have involved firstly, an intermediate mass black hole 

numbers of very bright X-ray sources external to the nu- binary. Makishima et al. (2000) required black hole masses 

clear regions, i.e. 36 super-Eddington sources of which 16 of between ~3 and ~80 Mq for the Eddington limit not 

were ULX (Fabbiano 1989). A significant fraction of the X- to be exceeded. In general, intermediate mass black hole 

ray emission from galaxies has long been known to originate models involve masses of 10 2 - 10 4 Mq, i.e. more massive 

from X-ray binaries (e.g. Fabbiano 1989; Fabbiano, Kim & than Galactic BHB such as CygX-1, and substantially less 

Trinchieri 1994, Fabbiano et al. 2001, Blanton et al. 2001). In massive than AGN. This possibility has been invoked to ex- 

a survey of nearby galaxies with Rosat, Roberts & Warwick plain the most luminous, variable ULX in M82 (Kaaret et al. 

(2000) found 28 ULX outside the nuclei. Similarly, Colbert 2001). However, King et al. (2001) discussed the difficulties 

& Mushotzky (1999) investigated extra-nuclear (>2 arcmin of forming these, and proposed an alternative model involv- 

offset) X-ray sources in 39 nearby galaxies finding 13 ULX. ing mild beaming and a link with Galactic micro-quasars. It 

They suggested these were accreting black hole systems of has been known that ULX occur preferentially in regions of 

10 2 -10 4 Mq. The spectra of 7 ULX in nearby galaxies stud- star formation (Zezas et al. 1999; Roberts & Warwick 2000; 

ied using ASCA by Makishima et al. (2000) were well-fitted Fabbiano et al. 2001), and the model is consistent with the 
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Table 1. Galaxy sample and Chandra observation log. ACIS exposures are shown, and the min- 
imum resolvable size scale (equivalent to 0.5 arcsec) using our best distances for the sources 
(Table 2). We also show the line-of-sight column densities within our galaxy for the centroid of 
each external galaxy. (Dickey & Lockman 1990). 



Galaxy 


Type 


T 




D 25 


jyGAL 


ObsID 


Date 


Exp 


Resin 










I 


10 20 cm" 2 




dd/mm/yr 


ks 


pc 


NGC4636 


E0-1 


-4.8 


10. 


6.0 


1.8 


323 


26/01/00 


52 


36 


NGC1132 


E 


-4.8 


13. 


2.5 


5.2 


801 


10/12/99 


13 


170 


NGC4697 


E6 


-4.7 


10. 


7.2 


2.2 


784 


15/01/00 


39 


28 


NGC 1399 


El 


-4.5 


10. 


6.9 


1.3 


319 


18/01/00 


56 


50 


NGC1291 


SO/a 


0.1 


9.3 


9.8 


2.1 


2059 


07/11/00 


23 


30 


NGC 2681 


SO/a 


0.4 


11. 


3.6 


2.4 


2060 


30/01/01 


77 


42 


NGC 253 


Sc 


5.1 


7.1 


28 


1.4 


969 


16/12/99 


14 


7.5 


NGC 3184 


Scd 


5.9 


10. 


7.4 


1.1 


804 


08/01/00 


41 


17 


NGC 4631 


Sd 


6.5 


8.6 


15 


1.3 


797 


16/04/00 


59 


8.5 


IC 5332 


Sd 


6.8 


11. 


7.8 


1.4 


2066 


02/05/01 


52 


9.7 


IC 2574 


Sm 


8.9 


10. 


13 


2.4 


792 


07/01/00 


8.6 


8.8 


NGC 1569 


Im 


9.6 


9.4 


3.6 


2.2 


782 


11/04/00 


92 


4.1 


IZW18 








0.30 


1.9 


805 


08/02/00 


31 


31 



expected association of high mass X-ray binaries (HMXB) 
with young stellar populations. The numbers of ULX found 
in the Antennae galaxies (Fabbiano et al. 2001) supports 
the connection with recent massive star formation. ULX 
were also detected in elliptical galaxies which do not contain 
HMXB, e.g. Sarazin et al. (2001), which led King (2002) to 
extend the model by proposing two types of ULX: persistent 
sources predominating in galaxies with young stellar popu- 
lations, and micro-quasars with bright, prolonged outbursts 
occurring in elliptical galaxies. 

However, there is also the possibility that ULX do not 
consist of a single or even two types of object. Objects gener- 
ating X-ray luminosities larger than the Eddington limit for 
a neutron star include not only stellar mass black holes (M 
< 100 Mq) but very young supernova remnants (Roberts & 
Warwick, 2000; King et al., 2001). Some pre- Chandra ULX 
may have been unresolved "super-bubbles" of shock-heated 
HII in the ISM (e.g. Stewart & Walter 2000), typically with 
diameters ~200-1000 pc. With the 0.3 arcsec (on-axis) res- 
olution of Chandra (van Speybroeck et al. 1997; van Spey- 
broeck 1999) these objects would appear extended in nearby 
galaxies, for example, 0.3 arcsec corresponds to a size of 27- 
40 pc at the distance of the Antennae galaxies (20-30 Mpc) 
(Fabbiano et al. 2001). Young, compact supernova remnants, 
which can reach luminosities of a few x 10 40 erg s _1 (Immler 
& Lewin 2002) would not be resolved even with Chandra, 
except in the nearest members of the Local Group. Simi- 
larly an unresolved cluster of sub-Eddington low mass or 
high mass XRB may appear ultra-luminous. This could be 
in a globular cluster such as the two sources in the Galactic 
globular cluster M15, White & Angelini 2001). Since the di- 
ameter of a globular cluster is typically 10 pc, these objects 
also cannot be resolved, even with Chandra. Large ampli- 
tude variability would be a strong pointer that the source 
consists of a single object or a group of a small number of 
objects only (Roberts & Warwick, 2000). 

In this paper, we examine systematically Chandra ACIS 
observations of 13 normal galaxies to detect all super- 
Eddington sources and so determine whether the prevalence 



of ULX depends upon galaxy morphology. One aim of the 
work was to examine the spectra of the brighter sources and 
search for any spectral property that may reveal the nature 
of a ULX source. 



2 THE GALAXY SAMPLE 

Thirteen galaxies were selected from the available Chandra 
observations in the public archive in the summer of 2001, 
that were classified by the Chandra X-ray Centre ( CXC) as 
Normal Galaxies. These galaxies were selected on the ba- 
sis that the observations were longer than ~10 ksec, and 
covered a range of types. While this is not a statistically 
complete sample, the seven spiral, four elliptical and two ir- 
regular galaxies allow us to investigate whether there is any 
correlation of ULX occurrence with morphology. The galax- 
ies chosen are listed in in morphological order from E to Irr 
in Table 1, with details of the observations. Morphological 
types and face-on, absorption-corrected B-band magnitude, 
Bj are taken from the Third Reference Catalogue of Bright 
Galaxies (RC3; de Vaucouleurs et al. 1991). The morpho- 
logical parameter T, where -5 is Elliptical, is SOa and 10 
is Irregular, is taken from the Lyon-Meudon Extragalactic 
Database LED A (Paturel et al. 1997), based on RC3. 

Chandra point-source detections have already been re- 
ported for NGC 1399 (Angelini, Loewenstein & Mushotzky 
2001), NGC 4697 (Sarazin, Irwin & Bregman 2000, 2001) 
and the central region of NGC 1291 (Irwin, Sarazin & Breg- 
man 2001). However, since spectral studies of individual ob- 
jects were not in general made, these sources were included 
here. 

Based on a literature search, we obtained preferred dis- 
tances for the 13 galaxies as shown in Table 2. In the case of 

5 galaxies, we use the distances given by Shapley, Fabbiano 

6 Eskridge (2001). For the other galaxies, we adopt, in or- 
der of preference the values obtained from Cepheid variables; 
from /-band surface-brightness fluctuations; from the Tully- 
Fisher relation; or from the isophotal diameter - velocity 
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Table 2. Preferred values of distances of each galaxy in Mpc 



Galaxv 


distan 


■ 

trCciinic|uc 


reference 


NGC 4636 


14.7 ±0.9 


SBF 


1 


NGC1132 


69 ±5 


D N -cr 


2 


NGC 4697 


11.7 ±0.8 


SBF 


1 


NGC 1399 


20 ± 1 


SBF 


1 


NGC 1291 


12 ±3 


D N -a 


3 


NGC 2681 


17 ±3 


SBF 


1 


NGC 253 


3.1 ±0.7 


Ccphcid 


4 


NGC 3184 


7.2 ± 1.7 


T-F 


4 


NGC 4631 


3.5 ±0.3 


T-F 


4 


IC 5332 


4.0 ±1.0 


T-F 


5 


IC 2574 


3.6 ±0.3 


Ccphcid 


4 


NGC 1569 


1.7 


BS 


4 


IZW18 


12±2 


v vir 


6 



The techniques shown are: SBF: /-band surface brightness 
fluctuation; D^-a: isophotal diameter velocity dispersion; T- 
F: Tully-Fishcr; BS: brightest stars in galaxies; V v i r : rcdshift 
distances with Virgoccntric correction. 

References: (1) Tonry et al. 2001; (2) Djorgovski & Davis 1987; 
(3) Prugniel & Simien 1996; (4) Shaplcy ct al. 2001; (5) Bot- 
tinclli 1985; (6) Ostlin 2000; 

dispersion (D^-a) relation. Otherwise, we adopt estimates 
based on the brightest star in the galaxy, or use the red- 
shift corrected for Virgocentric flow from LED A. Distance 
errors are also shown, reflecting the uncertainty in the most 
reliable technique available for any galaxy. Where possible, 
systematic and statistical errors are combined. If no errors 
were quoted in the literature we adopt, if possible, typical 
values for the given distance method as listed in Jacoby et 
al. (1992). For the redshift-distances we adopt errors deter- 
mined by Shapley et al. (2001) of -30% (20% in the case of 
NGC 1132). 



3 X-RAY SOURCE DETECTION 

Data analysis was performed using the CIAO 2.1.2 software, 
XANADU and FTOOLS 5.0. To remove periods of high back- 
ground, lightcurves were accumulated from source-free re- 
gions of the active chips and intervals having factor of two 
increases in count rate removed. This led to significant data 
loss only for IZW 18 and NGC 1291, in which ~25% and 
~35% of data were removed, respectively. Net exposure 
times are shown in Table 1. 

Searches for point sources were made using the CIAO 
wavdetect algorithm, without regard to whether detec- 
tions were at the centres of galaxies or not. However, as 
discussed below, only a few detections coincided with galaxy 
centres. The spurious source detection probability was set at 
10~ 6 per pixel, corresponding to 0.014 spurious detections 
per square arcmin. Exposure maps created with mkexpmap 
at 1.7 keV were used with wavdetect to minimize spurious 
detections at the detector chip boundaries and these aspects 
were checked by manual inspection of the images. In two 
galaxies, IC 5332 and NGC 1569, no sources were detected 
exceeding 5 x 10 38 erg s _1 . Detections spanning 2, 4, 8 or 
16 pixels were flagged, as these imply extended sources. Ra- 
dial profiles (Sect. 3.1) confirmed possible extended sources 



as those significantly broader than the PSF. For a typical 
off-axis PSF of 1 arcsec width, a typical distance of 10 Mpc 
implies a size of 50 pc, so that to appear point-like in Chan- 
dra a source must have a size no larger than this. Only 10 
sources eventually proved to have an extended component. 

We selected the point sources lying within the fJ-band 
25th magnitude isophotes derived from RC3, which are thus 
likely to be associated with the galaxy. For these, spectra, 
lightcurves and instrument response functions were gener- 
ated. Background data were extracted from annuli centred 
upon the source region and lying entirely within the same 
CCD node. The annuli were chosen to cover an area at least 8 
times larger than the source extraction region and contain- 
ing at least 20 photons (although such regions frequently 
exceeded this limit). To prevent contamination of the back- 
ground spectra by photons from bright point sources, all 
photons detected within a region six times larger than the 
1-a encircled-energy ellipse of every source, and centred on 
that source were removed prior to background accumula- 
tion. Such a large masking region was adopted since with a 
3-cr region, a few per cent of a source's intensity will remain. 

The total count from a given detection varied between 
~100 and 1000 (Table 4). Thus to plot light curves in the 
normal way to achieve Poisson errors of 10%, say, would al- 
low a very small number of bins, and searching for variability 
by x 2 testing would not be sensible in most cases. We thus 
applied a Kolmogorov-Smirnov test to the arrival times of 
photons and the numbers of photons that arrive by a par- 
ticular time are compared with those expected for constant 
intensity. This test strictly requires unbinned data, and so 
we used the primitive binning of the ACIS data of 3.24 s, 
although binning does not bias the test provided there are 
many primitive bins. This revealed only one source in galaxy 
NGC 253 which displayed variability at the 99% significance 
level, this source having a luminosity of 2.5 x 10 38 erg s _1 , 
located at a = 00 h 47 m 30. s 9, 8 = -25° 18' 26" (2000). Other 
sources exhibited no evidence for variability, however the 
Kolmogorov-Smirnov test does not provide upper limits for 
non-detections. We carried out simple \ 2 testing which in- 
dicated that other sources were not variable at a confidence 
of at least 90%. We adopt the procedure (below) of consid- 
ering only sources brighter than 5 x 10 38 erg s _1 , so that no 
bright sources detected (shown in Tables 3-5) displayed 
significant variability. Thus, none of the sources in the sam- 
ple examined provided evidence in this way that the object 
was not a superposition of sources. 

Radial brightness profiles were derived for each point 
source detected, which allowed extended sources to be 
clearly identified by comparing the profiles with the PSF 
derived for each position (see Sect. 3.1). Ten sources found 
using wavdetect were rejected as being clearly extended, 
and 2 more sources found to have point source and extended 
components are discussed in Sect. 3.1. 

Spectral fitting (Sect. 3.3) provided best-fit models from 
which fluxes were derived. Then, using the preferred dis- 
tances of Table 2, luminosities were calculated in the band 
0.3-7.0 keV. In the case of the weaker sources, which do 
not permit sensible fitting to discriminate between models, 
an absorbed power law was fitted with a fixed power law 
photon index of 2.0 (typical of fit results in other cases) 
and column density fixed at the line-of-sight Galactic value 
(Dickey & Lockman 1990; see Table 1). After removal of 



4 P. J. Humphrey et al. 



extended sources, in this band, 258 sources were found 
with L x > 1 x 10 38 erg s" 1 (Li), of which 158 (N 2 ) had 
L x > 2 x 10 38 erg s" 1 (L 2 ), 61 (N 5 ) had L x > 5 x 10 38 
ergs" 1 (L 6 ), and 22 (N m ) had L x > 1 x 10 39 erg s" 1 (L m ). 

In Table 3 we show the number of sources detected in 
each galaxy, N 2 , N 5 and N 10 . N 2 varied from 2 in NGC 1132 
to 72 in NGC 1399. The numbers detected have error values 
calculated using the errors in luminosity from spectral fitting 
combined with the distance uncertainties. In Fig. 1, Digi- 
tal Sky Survey images of each galaxy are shown with 25 th 
magnitude (D25) B-band isophotes from RC3 superimposed. 
All sources with luminosity > L w are shown, together with 
sources lying between L5 and L10. 

To estimate how many background objects may have 
been detected as ULX, the LogN(>S) - LogS relations in 
the Chandra Deep Field South (Giacconi et al. 2001) in the 
energy bands 0.3-2.0 keV and 2.0-7.0 keV were used to ob- 
tain the expected number of background objects within each 
galaxy (i.e. within the 25th magnitude isophote). Data in 
these bands were combined into a single plot for the total 
band, approximately, by using the spectral form of Giac- 
coni et al. (2001) to relate the total flux in the 0.3-7.0 keV 
band to the fluxes in the sub-bands. It was assumed that 
soft sources were also detected in the hard band, and so the 
correction was made using N values from the hard band and 
correcting S values to the flux in the total band. The number 
of background objects detected within each galaxy is shown 
in Table 3. For the 13 galaxies, 17 spurious sources are ex- 
pected with apparent luminosities > 5 x 10 38 erg s _1 , and 8 
with L x > lxlO 39 ergs" 1 . Hence 62±7% of the ULX sources 
are true members of their host galaxies (Table 3). To esti- 
mate the errors of the numbers of sources in each luminosity- 
band, the errors in distance and flux must be combined, 
bearing in mind that the distance errors do not affect each 
source independently. A simple Monte-Carlo method was 
adopted with 10000 simulations for each galaxy. Each sim- 
ulation consisted first of taking a random distance to the 
galaxy, normally distributed about the measured value, with 
a standard deviation equal to the 1-a distance uncertain- 
ties in Table 2. From this distance, appropriate flux limits 
for each luminosity-band were computed. A similar proce- 
dure was adopted taking a random flux for each detection 
distributed about its measured value, and the number of 
sources exceeding each flux limit were counted. We adopted 
the standard deviation of these numbers, accumulated over 
all the simulations, as the error estimate. 

We show the ULX detected in Fig. 1 (as diagonal 
crosses) in each galaxy, and other super-Eddington sources 
(L x > Lb) are shown separately. In Table 4 (upper panel), 
individual ULX in each galaxy are listed, and given a source 
name such as NGC 4636 PSX-1, the "P" indicating pointlike 
nature. For each source, the total number of counts ('count') 
contained in the spectrum (discussed in Sect. 3.3), the source 
count-rate, the luminosity (L x ) in the band 0.3-7.0 keV de- 
rived from spectral fitting using the best distance value, 
with 90% confidence errors, the offset (AR) of the source 
from the nominal galaxy centre. Source names following the 
CTiandra-naming convention are given. In the upper panel 
we include sources marked "-" which would fall below L10 at 
their lower limit. The lower panel shows sources exceeding 
Lb which would join the ULX detections if given their upper 



error limit luminosity (marked "+"). Table 5 similarly lists 
the extended sources detected. 

Next, we estimate the completeness of our source sam- 
ples, for each galaxy, i.e. the probability that a source of 
given luminosity will be detected by the detection algo- 
rithm. This requires firstly the estimated number of counts 
of the source in the image, this requiring a flux to counts 
conversion factor. We chose a conservatively low factor, i.e. 
an estimate erring on the low side of the counts obtained, 
by adopting a simple absorbed power law model having a 
photon index fixed at 1.0, i.e. a smaller value than typical 
of spectral fitting results. The column density was fixed at 
the average value measured in the sources detected. Aver- 
aging the threshold of the "correlation parameter" used by 
wavdetect to identify sources in the vicinity of each of 
our detections and assuming a Gaussian distribution of the 
correlation parameter (Freeman et al. 2002), we can esti- 
mate the fraction of the complete number of sources above 
each limit that will be detected by the algorithm. This is 
>97% for Lb sources (except in NGC 1132), and > 94% for 
L 2 sources (except for NGC 1132, NGC 1399, NGC 1291 and 
IZW 18). In these objects, the number detected represents a 
lower limit. 

As the Eddington limit for a 1.4 Mq neutron star is be- 
tween 2 — 4 x 10 38 erg s _1 , depending on composition, opac- 
ity and gravitational redshift (Paczyhski 1983), luminosities 
smaller than L5 do not exceed the Eddington limit substan- 
tially, and may be neutron star binaries. Only a very small 
number of Galactic neutron star binaries (< 1%) are known 
to ever exceed the Eddington limit: the bright LMXB GX 5-1 
has been observed with L x ~ 4 — 5 x 10 38 erg s _1 (Church & 
Balucinska-Church 2001; Christian & Swank 1997). In Sco 
X-l, the total luminosity often exceeds Lio- In Table 3 we 
show the detections of objects brighter than various thresh- 
olds. However, in the rest of the paper we concentrate on 
the ULX (L x > L10), and do not consider at all objects 
fainter than L5 thus excluding most neutron star binaries. 
The number of ULX detected is 22, and these exceed the 
Eddington limit for a 10 Mq black hole. There are 39 other 
sources in the luminosity range 5 — 10 x 10 38 erg s _1 which 
may be expected to consist mostly of Black Hole Binaries 
(BHB). 

3.1 Radial intensity profiles 

The excellent spatial resolution of Chandra (Weisskopf et al. 
2002) makes it possible to test for spatial extension of the 
ULX. Accordingly, radial intensity profiles for each detected 
bright source were extracted, to search for excesses above the 
point-spread function. Since the PSF width depends on the 
off-axis angle and the source spectrum, PSFs were generated 
for each source based on the mean detected photon energy 
using the CIAO tool MKPSF. Source data were extracted in 
18 logarithmically-spaced radial bins in the total band 0.3- 
7.0 keV and rebinned to ensure a minimum of 20 counts per 
bin, and compared with the PSF. This was done using dedi- 
cated software in which it was tested whether a point source 
model fits the data well. It was also possible to fit extended 
source models such as a simple King model. The background 
was assumed to be constant over the radial distance in the 
radial profiles, and all point sources detected in the neigh- 
bourhood of the source being modelled were excluded from 
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NGC 4636 EO-1 NGC 1132 E NGC 4697 E6 NGC 1399 El 




NGC 1291 SO/a NGC 2681 SO/a NGC 253 Sc NGC 3184 Scd 




NGC 4631 Sd IC 5332 Sd IC 2574 Sm NGC 1569 Im 




IZW 18 



Figure 1. Digital Sky Survey images of each of the galaxies in the sample with X-ray source detections within the £?-band 25th magnitude 
isophotes superimposed. In all cases, the grid unit is 1.2' square, except for IZW 18 which has a grid unit of 0.24' in Right Ascension 
and 0.18' in Declination, and the image has been magnified by a factor of 8 for clarity. The ULX (with L x > Lio) are indicated with an 
"X"; other sources (L x > L5) exceeding the Eddington limit for a neutron star are shown by a "+". Elliptical galaxies are shown at the 
top, then spiral galaxies, and finally the irregular galaxy. Except for NGC 253 and IC 2574, the detector chips enclosed all of the region 
within the isophote; in the latter galaxy the chip boundaries are shown. In IC 5332 and NGC 1569, no detections were made. 
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Table 3. The numbers of bright sources detected (see text); JV10 are the numbers of ULX. The number of sources with luminosities 
greater than 2 X 10 38 erg s" 1 (iV 2 ), 5 X 10 38 erg s _1 (A%), 10 39 erg s _1 (Nio) and 10 40 erg s _1 (iVioo) arc shown separately. Errors take 
into account both flux and distance uncertainties. The expected number of background sources (N bk9 ) for each lower limit luminosity is 
also shown, computed at the best distance. 



Galaxy 


Dist (Mpc) 


N 2 


JV 2 bkg 


N 5 


^bkg 


N 10 




iVioo 


IV 100 


NGC 4636 


14.7 ±0.9 


20 ±5 


3.7 


3 ±0.3 


1.5 


< 0.2 


0.74 





0.074 


NGC1132 


69 ± 5 


2 


9.9 


2 


4.0 


2 


2.0 


1 


0.2 


NGC 4697 


11.7 ±0.8 


14 ±2 


2.8 


4± 1.2 


1.1 


< 0.06 


0.57 





0.06 


NGC 1399 


20 ±1 


74 ±5 


11. 


26 ±3 


4.5 


9 ±2.5 


2.3 





0.23 


NGC 1291 


12 ±3 


21 ± 10 


7.3 


8 ±4.6 


3.0 


3± 1.8 


1.5 


1 


0.15 


NGC 2681 


17 ±3 


9 ± 1.2 


2.2 


5± 1.5 


0.87 


3 ±0.9 


0.43 





0.04 


NGC 253 


3.1 ±0.7 


10 ±3 


1.1 


7 ±2.2 


0.45 


3 ±2.0 


0.22 





0.22 


NGC 3184 


7.2 ±1.7 


3 ±2.8 


1.6 


2± 1.1 


0.65 


< 0.5 


0.33 





0.03 


NGC 4631 


3.5 ±0.3 


2 ±0.4 


0.31 


2 ±0.2 


0.12 


1±0.3 


0.062 





6 x 10~ 3 


IC 5332 


4.0 ±1.0 


< 0.5 


0.47 


< 0.2 


0.19 


< 0.06 


0.094 





9 x 10~ 3 


IC 2574 


3.6 ±0.3 


2 ±0.5 


0.50 


1 ±0.2 


0.20 


< 0.06 


0.11 





0.011 


NGC 1569 


1.7 ±0.2 


< 0.1 


9.7 x 10~ 3 





4.2 x 10" 3 





2.3 x 10" 3 





2 x 10~ 4 


IZW18 


13 ± 2 


1 ± 0.01 


6.0 x 10~ 3 


1 ±0.06 


2.4 x 10" 3 


1 ±0.3 


1.2 x 10~ 3 





io~ 4 


Total 




158 ± 13 


40.9 


61 ±6 


16.6 


22 ±4 


8.40 


2 


1.0 



Chandra ACIS 



NGC 253 PSX-3 



Chandra ACIS 



NGC 1291 PSX-1 




1 2 5 

Angular distance (arcsec) 



Angular distance (arcsec) 



Figure 2. Point sources embedded in extended emission. Radial brightness profiles of NGC 253 PSX-3 (left panel) and NGC 1291 PSX-1 
(right panel), showing the clearly extended emission. This has been parameterized in both cases using King models. The contributions 
of individual components in the fitting (background, PSF and King model) are shown. 



the background data. In the majority of sources, there was 
no evidence of extended nature. The extended sources are 
not included in Table 4, or in discussion of ULX, but arc 
shown separately in Table 5. 

In the case of the detection NGC 253 PSX-3, the ra- 
dial profile shown in Fig. 2 (left panel) consists of clear ex- 
tended emission many times wider than the PSF, together 
with excess intensity at the centre with a width comparable 
to the PSF, suggesting a point-source, and so this source 
in included in Table 4. In NGC 253 and NGC 4636, there 
are known to be regions of substantial, diffuse galaxy emis- 
sion exhibiting complicated structure (Strickland et al. 2000; 
Jones et al. 2002), so that it is difficult to test if sources close 
to the centres are real or concentrations of structured emis- 
sion. Similarly, where this extended emission has an erratic 



radial profile that is not easily fitted, it is difficult to de- 
termine whether sources superimposed upon it are broader 
than the PSF or genuine point-sources. In such ambiguous 
cases the objects are included in Table 5 and assumed to be 
extended. 

In five galaxies, detections were made within a few arc- 
seconds of the galaxy centres where there may be diffuse 
emission, so that the sources may not be point-like. The 
radial profiles of four of these sources (NGC 1132 ESX-1, 
NGC 1399 ESX-1, NGC 2681 ESX-1 and NGC 4697 ESX-1) 
could be fitted by King models, and so the sources are shown 
in Table 5 of extended sources, also containing extended 
sources offset from the centres as revealed by their radial 
profiles. The fitting results allowed the extended emission to 
be subtracted in fitting the radial profiles of other sources in 
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Table 4. Upper tabic: the ULX detections (including sources in which the errors would allow the luminosity to fall just below Lio, 
marked "-"); lower table: sources with luminosity close to Lio in which the luminosity errors would allow to join the ULX category (+ 
cases). The Chandra name includes the right ascension and declination of each source (epoch 2000); the total count, count rate and offset 
from the centre of the galaxy in arcseconds are shown. 



Source 


name 




count 


rate 
10~ 3 c s- 1 


L(0.3 - 7.0) 
10 38 erg s" 1 


AR 

a 


± 


ULX 


NGC 1132 
















PSX-1 


CXOU 


J025251.4-011631 


21 


1.6 


- .^,,-1-70 

1201™ 


4 




PSX-2 


CXOU 


J025249.4-011620 


10 


0.76 




36 




NGC 1399 














PSX-1 


CXOU 


J033832.6-352705 


700 


13.0 


41. ±4. 


54. 




PSX-2 


CXOU 


J033831. 8-352604 


570 


10.0 


22. ± 2. 


68. 




PSX-3 


CXOU 


J033820. 1-352446 


410 


7.3 


15. ± 2. 


190 




PSX-4 


CXOU 


J033827.6-352648 


520 


9.3 


14. ± 2. 


22 




PSX-5 


CXOU 


J033821. 9-352928 


190 


3.4 


13. ± 2. 


180 




1 OVV-U 




JUOOozy. 1 - JOjZOWi 


9^n 
zou 


4 1 


lO. Q 


1 1U 




PSX-7 


CXOU 


J033825.2-352522 


340 


6.1 


11. ± 1. 


110 




PSX-8 


CXOU 


J033832.3-352710 


120 


2.1 


11 + 5 ' 


51 




PSX-9 


CXOU 


J033833. 1-352731 


230 


4.1 


11. ± 2. 


70 




NGC 1291 
















PSX-1 


CXOU 


J031718.6-410629 


390 


17.0 


130+140 


4 




PSX-2 


CXOU 


J031702.5-410714 


360 


16.0 


18. ± 2. 


230 




PSX-3 


CXOU 


J031713.8-411035 


310 


13.0 


17. ± 3. 


250 




NGC 2681 
















PSX-1 


CXOU 


J085335.7+511917 


210 


2.7 


44+|; 


46 




PSX-2 


CXOU 


J085324.4+511819 


45 


0.6 


18+| 


130 




PSX-3 


CXOU 


J085333.7+511930 


380 


4.9 


12. ± 1. 


37 




NGC 253 
















PSX-1 


CXOU 


J004734.0-251637 


650 


46.0 


12. ± 1. 


44 




PSX-2 


CXOU 


J004733.0-251749 


1100 


79.0 


12.1 ±0.8 


31 




PSX-3 


CXOU 


J004733.4-251722 


810 


58.0 


11 + 5 ' 


6 




NGC 4631 
















PSX-1 


CXOU 


J124211. 1+323236 


1100 


19.0 


13.4 ±0.9 


52 




IZW18 
















PSX-1 


CXOU 


J093402.0+551428 


390 


13.0 


16. ±2. 


5 




Possible ULX detections 


NGC 4636 
















PSX-1 


CXOU 


J124249.1+024046 


140 


2.7 




33 


+ 


PSX-3 


CXOU 


J124255.8+024035 


69 


1.3 




99 


+ 


NGC 1399 
















PSX-10 


CXOU 


J033828.6-352724 


58 


1.0 




27 


+ 


PSX-11 


CXOU 


J033826.5-352732 


150 


2.7 


9.5 ±2.1 


50. 


± 


PSX-12 


CXOU 


J033832.8-352658 


140 


2.5 


9.2 ±2.9 


56. 


± 


PSX-13 


CXOU 


J033827.9-352747 


78 


1.4 


9.0+H 


51 


+ 


PSX-14 


CXOU 


J033832.3-352702 


98 


1.8 


8.8+J;* 


50 


+ 


PSX-15 


CXOU 


J033815.4-352628 


210 


3.7 


8.5 ±1.9 


200 


+ 


PSX-17 


CXOU 


J033833.8-352658 


71 


1.3 


7 3+11. 
'•°-4.1 


71 


± 


PSX-18 


CXOU 


J033827.8-352750 


110 


2.0 


7 q+16. 
' -°-2.4 


55 


± 


PSX-22 


CXOU 


J033830.4-352430 


90 


1.6 


6 4+ 9 ' 2 


150 


+ 


NGC 1291 














PSX-4 


CXOU 


J031712.1-410438 


27 


1.2 


7.8 ±5.3 


140 


+ 


NGC 2681 
















PSX-4 


CXOU 


J085336.4+511727 


59 


0.8 


8.2 ±2.9 


99 


+ 
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Table 4. contd. Other super- Eddington detections (L x > L5). In none of these sources do the luminosity errors 
allow the source to join the ULX category 



Source 


name 




count 


rate 
10~ 3 c s- 1 


L(0.3 - 7.0) 
10 38 erg s" 1 


AR 

// 


Other super-Eddington detections 


NGC 4636 














PSX-2 


CXOU 


J124251.4+023835 


220 


4.2 


7.0±H 


160 


NGC 4697 














PSX-1 


CXOU 


J124846.8-054854 


350 


9.0 


7.0 ±1.2 


170 


PSX-2 


CXOU 


J124833.2-054742 


200 


5.1 


6.4 ±1.1 


45 


PSX-3 


CXOU 


J124836.7-054732 


130 


3.3 


5.3 ± 1.2 


32 


PSX-4 


CXOU 


J124839.3-054808 


190 


4.9 


5.2 ±0.8 


51 


NGC 1399 














PSX-16 


CXOU 


J033831. 9-352649 


170 


3.0 


7.5 ± 1.1 


44. 


PSX-19 


CXOU 


J033831. 3-352411 


200 


3.6 


7.3 ±1.6 


170 


PSX-20 


CXOU 


J033828.9-352602 


130 


2.3 


7.1 ± 1.7 


56 


PSX-21 


CXOU 


J033836.8-352747 


200 


3.6 


7.0 ± 1.0 


130 


PSX-23 


CXOU 


J033836.3-352809 


61 


1.1 


6.4 ±2.2 


130 


PSX-24 


CXOU 


J033827.2-352601 


110 


2.0 


5.2 ± 1.4 


63 


PSX-25 


CXOU 


J033831. 6-352600 


160 


2.9 


5.1 ±0.9 


70 


PSX-26 


CXOU 


J033816.5-352745 


60 


1.1 


— 2.. -j 


190 


NGC 1291 














PSX-5 


CXOU 


J031658.7-410406 


67 


2.9 


6 4+ 2 ' 4 
"•^-1.9 


310 


PSX-6 


CXOU 


J031718.9-410628 


79 


3.4 


5.4 ± 1.7 


20 


PSX-7 


CXOU 


J031721. 7-410653 


83 


3.6 


5.1 ± 1.4 


65 


PSX-8 


CXOU 


J031719.3-410615 


86 


3.7 


5.0 ± 1.4 


29 


NGC 2681 














PSX-5 


CXOU 


J085333.5+511852 


140 


1.8 


5.1 ± 1.2 


6 


NGC 253 














PSX-4 


CXOU 


J004717.6-251811 


770 


55.0 


7.8 ±0.6 


240 


PSX-5 


CXOU 


J004722.6-252051 


880 


63.0 


6.5 ±0.5 


260 


PSX-6 


CXOU 


J004732.1-251721 


140 


10.0 


5.7 ± 1.5 


16 


PSX-7 


CXOU 


J004735.3-251512 


490 


35.0 


5.2 ±0.5 


130 


NGC 3184 














PSX-1 


CXOU 


J101823.0+412742 


460 


11.0 


7.1 ±0.8 


160 


PSX-2 


CXOU 


J101812.0+412421 


420 


10.0 


5.0 ±0.6 


100 


NGC 4631 














PSX-2 


CXOU 


J124155.6+323217 


3200 


54.0 


7.2 ± 0.3 


180 


IC 2574 














PSX-1 


CXOU 


J102843.1+682816 


440 


51.0 


6.6 ±0.8 


350 



these five galaxies, and when this was done, no other source 
within these galaxies proved to be extended. In the case of 
NGC 1291, a point-source PSX-1 was detected at the galaxy 
centre, superimposed on extended emission (see Fig. 2, right 
panel); this is discussed in Sect. 5.3. 

3.2 Hardness ratios 

As a simple method of comparing the sources, hardness ra- 
tios were obtained for all sources in Table 4 based on three 
energy bands: a soft band (S): 0.3-1.5 keV, a medium band 
(M): 1.5-3.0 keV and a hard band (H): 3.0-5.0 keV. Hard- 
ness ratios were defined as HR1 = (S-M) /(S+M) and HR2 = 
(M-H)/(M±H). The HR1 ratio is more sensitive to changes 
in column density, while the HR2 ratio is more sensitive to 
the power law index. With this definition, a harder source 
implies a more negative value of HR2. A colour-colour dia- 
gram (HR2 versus HR1) is shown for the detected sources 
in each galaxy in Fig. 3, plotting each source as a point with 



error bars derived from the Poisson errors of the intensities 
in each energy band. 

To investigate quantitatively the sensitivity of the ratios 
of these bands to spectral changes, ACIS data were simu- 
lated for a simple absorbed power law spectral model for 
a range of column densities Nn and power law indices V. 
Count rates in each of the energy bands were calculated and 
points on the (S/M, M/H) colour-colour diagram found, so 
giving a T/Nh grid (c.f. Kim et al. 1992). This grid is super- 
imposed on the colour-colour diagram for the sources found 
in each galaxy in Fig. 3. (The V scale is vertical, the Nh 
scale horizontal, with increasing values to the top and left.) 

Although the spectra of LMXB, HMXB or BHB cannot 
generally be described by this simple one-component model, 
an absorbed power law model can be fitted to X-ray binary 
spectra as a method for gauging the steepness of the decrease 
to high energies, even if the quality of fit may be poor. For 
example, CygX-1 (a BHB) in the High State has a very 
soft spectrum because of the strong disc blackbody compo- 
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Table 5. The extended source L5 detections 



Galaxy 


source 


a 


<5 


count 


rate 


L(0.3 - 7.0) 


AR 






h m s 


0/ // 




10~ 3 c s- 1 


10 38 erg s _1 


tt 


NGC4636 


















ESX-1 


12 42 50.2 


2 41 15 


620 


12.0 


9 o+49 
z -13 


5.8 




ESX-2 


12 42 49.8 


2 41 15 


900 


17.0 


181? 


1.7 




ESX-3 


12 42 49.7 


2 41 11 


1100 


21.0 


16+4 


6.3 




ESX-4 


12 42 50.2 


2 41 10 


400 


7.7 


7.5 ±2.9 


9.4 


NGC1132 


















ESX-1 


2 52 51.8 


-1 16 28 


100 


7.7 


140 ± 30 


5.5 


NGC 4697 


















ESX-1 


12 48 35.9 


-5 48 3 


180 


4.6 


4.5 ±0.9 


1.0 


NGC 1399 


















ESX-1 


3 38 29.0 


-35 27 1 


11000 


200.0 


760 -80 


3.4 


NGC 2681 
















ESX-1 


8 53 32.7 


51 18 49 


900 


12.0 


19 ±2 


6.8 


NGC 253 


















ESX-1 


47 33.3 


-25 17 17 


410 


29.0 




2.4 




ESX-2 


47 33.0 


-25 17 20 


600 


43.0 


10 ± 1 


2.1 



nent (Dotani et al. 1997). Although the power law index of 
the Comptonized emission varies from ~1.7 to 2.0 as the 
source moves from the Low State to the High state, fitting 
the one-component model to High State data will yield an 
apparent T of ~3.5, compared with ~1.9 in the Low State. 
We have carried out simulations based on previous analysis 
of Galactic Black Hole binaries and Low Mass X-ray Bina- 
ries, and results are shown in Fig. 4. The simulation based 
on the CygX-1. High State data of Baluciriska-Church et al. 
(1998) confirms a location on the colour-colour diagram at 
a grid value with T = 3.5, so that sources located at sim- 
ilar positions in Fig. 3 may be Black Hole Binaries in the 
High State. Power law indices between 2.5-3.5 may indicate 
black hole nature; however given the errors of points in Fig. 
3, it would be unsafe to associate black hole nature except 
with values >3.0. There were five such objects, including two 
ULX: NGC 1399 PSX-1 and PSX-2, plus NGC 1399 PSX-15 
(luminosity marginally below 10 39 erg s" 1 ), and PSX-16 and 
PSX-26, having L x > L5. Note that large values of T may 
also result from soft thermal spectra, for example those of 
SNR or supersoft sources. 

Colour-colour diagrams however cannot distinguish a 
BHB in the Low State from a LMXB, although most LMXB 
should have been excluded from the sample defined with 
L x > 5 x 10 38 erg s _1 . Moreover, it would not be possible to 
distinguish different classes of LMXB; for example, the flar- 
ing branch in Z-track sources is distinguished from the nor- 
mal branch by flaring taking place above 7 keV (Baluciriska- 
Church et al. 2001), but tests show that this behaviour would 
not cause any movement in the present colour-colour di- 
agram. However, the diagram does allow comparisons be- 
tween galaxies, and with galaxy morphology, ft should be 
noted that fitting the absorbed power law model to ACIS 
data was adequate in only half of the sources shown in Table 
6 for which spectral fitting was carried out. 

Differences in hardness ratio between galaxies can be 
seen. In the early-type galaxies, the column densities are 
low: <5 5 x 10 21 cm~ 2 , as expected since early type galaxies 
generally have a hot ISM. Higher column-densities are seen 



however, in NGC 4636 PSX-3 and NGC 2681 PSX-1, PSX-2 
and PSX-4, sources that are not near the nucleus. Examina- 
tion of the spectra of these sources reveals evidence of heavy 
absorption in each case. This may be intrinsic, for exam- 
ple if the sources are strongly-absorbed background AGN. 
We note that the expected numbers of background objects 
(Sect. 3) is consistent with the numbers of heavily-absorbed 
sources we detect in elliptical galaxies. Alternatively, they 
may be intrinsically-absorbed objects within the host galaxy. 
In the nearly edge-on spiral galaxy NGC 253 we see evidence 
of a spread of column-densities consistent with the presence 
of cold ISM in spirals, reaching a few xlO 22 cm~ 2 . Spec- 
tral fitting results (Table 6) confirm high values of column 
density in the sources in Fig. 3 with large Ah . 



3.3 Spectral analysis 

We restrict spectral fitting to sources with L x > 5 x 10 38 
erg s _1 and for a total number of counts in the spectrum of 
> 200, which is a minimum for useful spectral fitting. This 
left a sample of 22 sources consisting of 14 ULX plus 8 other 
super-Eddington sources (L x > Lb). Spectral fitting results 
are shown in Table 6 for the ULX (upper panel), and other 
sources (lower panel). 

Spectra were extracted from the elliptical regions pro- 
vided by the detection algorithm and background data ex- 
tracted from annuli centred on each source from events files 
from which all point sources had been excluded. Data below 
0.3 keV, where the calibration becomes uncertain, and above 
7.0 keV, where the particle background exceeds the source 
count rate, were ignored. Spectra were grouped to a mini- 
mum of 20 counts per bin to allow use of the \ 2 statistic. 
However, the full energy range could not be used in many 
cases because of insufficient counts. Although dependent on 
the slope of the spectrum, to have usable data above 5 keV, 
more than 500 counts are needed. Non-standard spectral 
models were written for use within xspec in which the in- 
tensity or normalization of the model equals the unabsorbed 
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Figure 3. X-ray colour-colour plot for sources with L x > 5 X 10 38 erg s — 1 in each galaxy, ordered by morphological parameter T. ULX 
sources are shown with darker lines; other super-Eddington sources are shown with fainter lines. Each panel shows a grid of power law 
index Y - Nn values obtained by simulating ACIS data for an absorbed power law spectral model (see text). T increases vertically with 
values between - 3.0 in steps of 0.5, and N H increases to the left with values 0.0, 0.25, 0.5, 1.0, 2.0, 4.0 and 8.0 xl0 22 cm- 2 . In IC 5332 
and NGC 1569 no super-Eddington sources were detected, so they are omitted. 



energy flux in the 0.3-7.0 keV band thus allowing flux errors 
to be obtained directly as the error of a single parameter. 
Simple absorbed power law (P), multi-colour disc blackbody 
(D) and Mekal (M) models were tested against the data, 
as were two-component combinations of these. The Mekal 
model would be appropriate to optically thin emission with 
lines from hot plasma, appropriate for some SNR. In no cases 
was a cut-off power law representation of Comptonization 



fitted given that the spectra did not extend beyond 7 keV 
and use of this model is likely to result in curvature at ~5 
keV actually due to thermal emission being incorrectly in- 
terpreted as a high energy cut-off. In Table 5, results shown 
in boldface indicate that the model highlighted produced an 
appreciable decrease in \ 2 compared with other models. 

Of the 14 ULX in Table 6, 5 could be fitted by a sim- 
ple absorbed power law, although in none of these cases was 



ULX in Normal Galaxies 11 



i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



Cyg X-l High State 




-1.5 -1 -0.5 0.5 1 1.5 

HH1 

Figure 4. Calibration of grid of X-ray colour-colour diagram us- 
ing Galactic Black Hole Binaries and Low Mass X-ray Binaries 
(see text). CygX-1 and GRO J1655-40 are BHB, the other sources 
are LMXB. Simulations were based on the work of Dotani et al. 
(1997), Balucinska-Church & Church (2000), Balucinaka-Church 
et al. (2001), Church & Balucihska-Church (2001), and Smale et 
al. (2002). Column density and power law index contours corre- 
spond to the same values of these parameters as in Fig. 3 



the power law photon index large enough (r > 3.0) to be 
taken as evidence for black hole nature (Sect. 3.2). In one 
source (NGC 1399 source PSX-3) a Mekal model was pre- 
ferred marginally, offering some evidence for SNR nature. 
In two sources (NGC 253 PSX-1 and NGC 4631 PSX-1) an 
absorbed blackbody model was preferred although it was not 
possible to discriminate between disc blackbody and simple 
blackbody. These may indicate a black hole in a high state. 

The other 7 sources required two-component models, 5 
being fitted by a blackbody plus power law which would sug- 
gest neutron star binary or black hole binary nature. How- 
ever, it was not possible to discriminate between the simple 
blackbody and disc blackbody spectral forms and so discrim- 
inate between a neutron star binary and a black hole binary. 
IZW18 PSX-1 was best fitted by a Mekal plus power law, 
suggesting possible SNR nature, as was NGC 1291 PSX-1, 
although in this case it was necessary for the power law to 
have extra absorption. 

Comparison was next made of spectral fitting and hard- 
ness ratio results for ULX. Two ULX were identified from 
the harness ratio diagrams with a spectral shape which when 
fitted by an absorbed power law model, would give T > 3.0. 
These were NGC 1399 PSX-1 and PSX-2. Spectral fitting re- 
sults for the second of these had an actual power law photon 
index of 2.5-3.5 as part of a two-component model, adding 
some support to the idea that it is a High State BHB. 

Of the 8 other (non-ULX) sources shown in the lower 
panel of Table 6, exceeding the Eddington limit for a neu- 
tron star, 4 can be fitted by an absorbed power law model, 
although none of these had a photon index > 3.0 suggesting 
black hole nature. The other 4 were fitted by two-component 
blackbody plus power law models consistent with neutron 
star binary or black hole binary nature. In one of these, 
NGC 253 PSX-7, there was evidence for a steep power law 
index (> 3.0) strongly suggestive of black hole binary nature. 
The luminosity of this source is 5.2 x 10 38 erg s _1 . Very few 



Galactic neutron star binaries have been observed with such 
luminosities, and as the luminosity is substantially less than 
the Eddington limit of a 10 Mq black hole, the object is 
likely to be a black hole binary similar in the high state. 
Fig. 5 illustrates higher quality data, showing the best fits 
to two sources. NGC 1399 PSX-1 (left panel) is fitted by an 
absorbed power law model, whereas the second source, 
NGC 253 PSX-4 is well-fitted by a two-component model 
consisting of either a simple blackbody or a disc blackbody 
plus a power law. 



4 SEARCH FOR CORRELATION OF ULX 
NUMBERS WITH STELLAR MASS AND 
STAR-FORMING MASS 

We have investigated whether correlations exist between the 
numbers of ULX (iVio) and the numbers of super-Eddington 
sources (N5) with the stellar masses of the galaxies, and with 
the mass of star-forming regions in the galaxies. We used for 
the first test the near-infrared //-band fluxes accepted as a 
good measure of stellar mass, and for the second test, we 
used the 60 /im fluxes. In both cases, data were obtained 
from the NASA Extragalactic Database (NED) which also 
provided angular diameters for the galaxies in the form 61 
and 62 each given in arcminutes. For each galaxy, we used 
the total //-band flux / in Janskys from the photometric 
measurements listed in the database. We adopted the pro- 
cedure of taking f / 61-62 as a measure of the surface bright- 
ness of the //-band luminosity of each galaxy. The number 
counts iVio and N& were converted to counts per unit area 
by dividing by the projected area of each galaxy d 61 ■ 62, 
where d is the galaxy distance. It is well-known that this 
procedure avoids possible biases due to varying sizes of the 
galaxies, //-band fluxes were available from the database for 
all of the galaxies in the sample, except IZW 18. Similarly, 
there was no 60 /im flux for IZW 18 and only upper limits 
for NGC 1399 and NGC 1132. These upper limits were not 
included in Fig. 7 in which we test for correlation of number 
counts with 60 fim flux. We use the N10 and N5 number 
counts listed in Table 3 which include 5 upper limits for iVio 
and 1 upper limit for N$; these upper limits are plotted in 
Figs. 6 and 7 but are excluded from the correlation analysis. 
In NGC 1569, no sources were detected. 

In Fig. 6 we show the variation of N10 and N5 per unit 
area with //-band surface brightness as an indicator of the 
stellar mass of each galaxy. Neither of these suggests a cor- 
relation with stellar mass. The modified Kendal's-r statistic 
demonstrated a 88.1% probability of no correlation between 
iVio and stellar mass, and a 78.8% probability of no corre- 
lation between N5 and stellar mass. In Fig. 7 we show the 
corresponding plots for 60 /im surface brightness. The N5 
counts reveal a strong correlation with the infrared surface 
brightness with a correlation coefficient of 0.71 and proba- 
bility of 98.7% of correlation. In the case of the N10 data, 
there are only 4 points that are not upper limits, but these 
also can be seen to be correlated with 60 /im flux, although 
the formal probability of 82.6% may be regarded as insuffi- 
ciently high to prove the correlation. 

It is clear that the numbers of ULX and of super- 
Eddington sources correlate with 60 pm flux suggesting a 
connection with star formation. Fig. 7 shows that the two 
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Table 6. Spectral fitting results: upper panel: ULX; lower panel: other super-Eddington sources. Only spectra with > 200 counts 
were fitted (see text). The spectral models used are: P = power law, B = blackbody, D = disc blackbody, M = Mekal, in all cases 
with cool absorber. The M+AP model has additional absorption for the power law of 2.2 ± 1.0 ■ 10 22 atom cm . Fitting results in 
boldface identify the models which were significantly preferred. Although the source NGC 1399 PSX-9 had more than 200 counts, 
fitting had less than 8 dof and for this reason, sensible fitting was not possible. 
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Figure 5. Best fits to the spectra of the ULX sources NGC 1399 PSX-1 and NGC 253 PSX-4, chosen to illustrate higher quality data 
from the sample. NGC 1399 PSX-1 was fitted by an absorbed power law model, and NGC 253 PSX-4 by an absorbed blackbody plus 
power law, labelled B and P. 



elliptical galaxies plotted NGC 4636 and NGC 4697 appar- 
ently follow the same trend as the spiral galaxies, but with 
generally smaller number counts per unit area and the small- 
est infrared surface brightnesses. The amount of star forma- 
tion in elliptical galaxies would be expected to be smaller 
than in spiral galaxies, assuming that the infrared emission 
arises from star formation and not from an active nucleus. 
A full investigation of these and other individual sources in 
our sample will require extensive future work; however, we 
note that NGC 4636 is a low- luminosity AGN with radio jets 
and is a LINER (Nagar et al. 2000). NGC 4697 is a very flat- 
tened galaxy which is thought may, in fact, be of SO type 
(e.g. Dejonghe et al. 1996). Of the spiral galaxies in Fig. 7, 
NGC 2681 is a LINER in a SO galaxy with a circumnuclear 
ring of star formation (Gonzalez Delgado et al. 1997). The 
galaxy NGC 253 with substantial numbers of bright sources 
is a starburst galaxy (see Engelbracht et al. 1997) . 



5 DISCUSSION 

In a survey of 13 normal galaxies we have detected 22 ULX 
(with L x > Lio); of these, 15 exceed Lio by less than a fac- 
tor of two while the remainder have higher luminosities, 5 of 
which are « 10 40 erg s -1 . Thus the fainter objects are not 
inconsistent with the Eddington limit for a 10 Mq black 
hole binary, whereas the brighter objects exceed this limit 
by factors of 30 - 100. We have found no evidence of correla- 
tion with stellar mass but substantial evidence for a correla- 
tion with star forming mass. Fig. 7 shows that more bright 
sources are found in spiral galaxies consistent with the larger 
star forming mass in these. 



5.1 Nature of the ULX 

In many cases, our results from hardness ratio studies do 
not discriminate between types of object. Similarly, because 
of poor statistics, spectral fitting does not in many cases re- 
veal the nature of the object. Consequently, we concentrate 



discussion on those cases where the results do help identify 
the type of object. 

The colour-colour diagrams revealed 2 objects consis- 
tent with High State black hole binary nature: NGC 1399 
PSX-1 and PSX-2. Spectral fitting results were obtained for 
14 ULX. None of these required a photon index > 3 suggest- 
ing black hole nature; however, two sources were well-fitted 
by a one-component blackbody model (NGC 253 PSX-1 and 
NGC 4631 PSX-1), and this suggests BHB nature. The lumi- 
nosities of these two sources were 1.2 x 10 39 and 1.3 x 10 39 
erg s _1 , i.e. not exceeding the Eddington limit for a 10 Mq 
black hole. Two-component models similarly did not provide 
evidence for black hole nature. 

The correlation between ULX and super-Eddington 
source numbers and 60 fim flux supports previous results 
that such sources occur preferentially in star forming regions 
(Zezas et al. 1999; Roberts & Warwick 2000; Fabbiano et al. 
2001). Moreover, our results also show detection of sources 
in elliptical galaxies, and interestingly, Fig. 7 shows that N5 
in the two elliptical galxies NGC 4636 and NGC 4697 appear 
to follow the same correlation as for the spiral galaxies, sug- 
esting that star formation could be important in these two 
galaxies also. 

5.2 Comparison with previous ULX number 
counts 

Next, we compare the numbers of detections in the present 
work with previous surveys. Roberts & Warwick (2000) re- 
ported results from the Rosat HRI for a sample of 486 B- 
band bright northern galaxies, and detected 142 non-nuclear 
sources. We have used their luminosity distribution dN/dL3g 
= (1.0 ± 0.2) Ljjg 1 ' 8 , with L38 normalized to a B-band lumi- 
nosity of 10 10 Lq to calculate expected numbers of bright 
objects in our sample. In NGC 1399 and NGC 253 ULX 
counts of 1.3 and 0.4 are expected, and 7 and 3±2 respec- 
tively, obtained. In these two galaxies, N5 values of 2.2 and 
0.7 are expected, and 21 and 6.5 obtained. It appears that 
there is an underestimation of both ULX and N5 sources 
in the Rosat survey by a factor of 5 - 10 compared with 
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Figure 6. ULX count per unit projected area (left panel) and supcr-Eddington source count per unit projects area (right panel) versus 
_ff-band surface brightness. Open squares show elliptical galaxies and closed circles show spiral galaxies. 
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Figure 7. Nio per unit area (left panel) and W5 per unit area (right panel) versus 60 )im surface brightness. Open squares show elliptical 
galaxies and closed circles show spiral galaxies. 



Chandra, indicating the improved sensitivity and angular 
resolution of Chandra. 

Grimm et al. (2003) using previous studies of High Mass 
X-ray Binaries in the Milky Way and Magellanic Clouds, and 
Chandra and ASCA observations of nearby starburst galax- 
ies, showed a linear relation between the star formation rate 
(SFR) as measured by far infrared flux and HMXB collec- 
tive luminosity. They proposed that the number or collective 
X-ray luminosity of HMXB could thus be used as a measure 
of star formation rate. As an example, we will compare with 
the galaxy in our sample having the highest 60 /xm surface 
brightness, the starburst galaxy NGC253 (Fig. 7). We can 
adopt a value of SFR/Ma yn , where Md yn is the dynamical 
mass of the galaxy of ~ 2.0 x 10~ n yr _1 from Grimm et al. 
for galaxies of similar type (Sc). Taking M dyn as 1.2 x 10 11 
Mq (calculated from the maximum rotation velocity from 
LED A), this gives a value of SFR ~ 2.4 M yr _1 which 



combined with the linear relation (Fig. 3 of Grimm et al.) 
gives a count N2 = 7. This compares well with our measured 
value of 9 (Table 3). 

In the case of the Antennae galaxies, our results can 
be used to predict the expected numbers of bright objects. 
Using approximate fits to the data in Fig. 7 plus the 60 
/jm flux of the Antennae galaxies, we estimate a value of 
N 5 of 39, and a value of Nio of 21.7 agreeing reasonably 
with the 29 super-Eddington and 17 ULX actually detected 
(Fabbiano et al. 2001). 



5.3 The most luminous sources 

Two X-ray sources have luminosities exceeding 10 40 erg s _1 : 
NGC 1132 PSX-1 and NGC 1291 PSX-1. These luminosities 
imply a mass of at least 100 Mq if emitting isotropically, 
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assuming they are Eddington limited, and it is important to 
consider these cases carefully. 

The former source has 21 counts, so its luminosity is 
poorly determined, and the error-bars are sufficiently large 
that the source luminosity may be as small as 70 x 10 38 erg 
s _1 . This is still a very luminous object, ft is situated away 
from the nucleus of its parent galaxy, and we expect 0.2 
background sources of similar apparent brightness (Sect. 3), 
so we cannot exclude the possibility that it is not associated 
with the galaxy. NGC 1291 PSX-1 has a high count of 390, 
and a lower limit luminosity of 60 x 10 38 erg s _1 . It is a 
point source embedded in extended emission at the centre 
of its galaxy, and may be similar to low luminosity AGN-like 
cores in many nearby galaxies (Ho et al. 2001). 

5.4 NGC 1399 and Globular Cluster X-ray 
Sources 

NGC 1399 is the massive central galaxy in the Fornax Clus- 
ter having an extensive globular cluster system (Kissler- 
Patig ct al. 1999). Individual globular clusters will not gener- 
ally be resolved and so multiple sources within them may be 
confused. Indeed Angelini et al. (2001) note a higher average 
count rate of apparently point-like sources in globular clus- 
ters in NGC 1399, which they suggest may be explained in 
this manner. This may be the reason why there are so many 
ULX in this galaxy. Bright globular cluster X-ray sources are 
usually LMXB. Although there is some evidence for globular 
cluster sources being somewhat less bright on average than 
LMXB in the rest of the Galaxy, Verbunt (2002) concludes 
that this may be unreal. Thus if we use an average luminos- 
ity of non-cluster sources of 2 x 10 37 erg s _1 (Verbunt et al. 
1984), we find that the addition of 50 of these would needed 
to explain any ULX in this galaxy. 

5.5 Limits on Non- X-ray Binary Populations 

The super-Eddington sources and the ULX in our sample are 
point-like, within the spatial resolution of Chandra, implying 
sizes typically < 50 pc. These limits would seem to preclude 
super-bubbles or super-giant shells being possible sites for 
the emission. Although young SNR may be compact, we 
would expect these to be found preferentially in regions of 
continuing star formation in late-type galaxies. SNR spectra 
may either appear as optically-thin emission with lines, or 
alternatively be featureless. Our spectral fitting shows Mekal 
models as preferred in only one or two cases, suggesting a 
SNR nature; we may however detect other SNR with fea- 
tureless spectra which are not revealed in this way. Zezas 
ct al. (2002) show by optical identification that a relatively 
small fraction of ULX in their sample are SNR, implying 
that our dependence on infrared flux must be mostly due to 
other types of object. 



6 CONCLUSIONS 

This survey of 13 galaxies revealed 22 ULX sources with 
L x (0.3-7.0 keV)> 1 x 10 39 erg s~\ and a further 39 sources 
with L x > 5 x 10 38 erg s 1 (but less than 1 x 10 39 erg s _1 ). Ac- 
cidental coalignment of background sources is not a major ef- 
fect as 62% of these ULX must truly lie within their apparent 



host galaxies. Several sources were found that were clearly 
extended beyond the instrumental PSF and are shown in Ta- 
ble 5. Two sources appeared on the colour-colour diagram 
in a position suggesting possible black hole nature. 

We present spectral fitting results for 14 of the 22 ULX 
having a count > 200, the total count varying between 
200 and 3200 in the spectra fitted. Several of the sources 
could be fitted by an absorbed power law model suggest- 
ing non-thermal emission such as Comptonization. Seven of 
the sources were well-fitted by a blackbody plus power law 
model, also suggesting X-ray binary nature, although it was 
not possible to distinguish the type of XRB. Two sources 
were well-fitted by a one-component blackbody model sug- 
gesting High State black hole nature. However, we do not 
have positive evidence that any more than a relatively small 
fraction of the ULX are black hole binaries. Moreover, we 
find evidence that the numbers of ULX and super-Eddington 
sources are correlated with the star formation rate as indi- 
cated by the 60 fim flux supporting similar results found pre- 
viously, which is consistent with a link with younger stellar 
populations such as high mass X-ray binaries. 
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